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Abstract. In this paper, we develop a method based on the analysis of the 
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1.1. Homogeneous spaces 

The classification of liomogeneous and anisotropic 3D spaces gives 11 different types of 
spaces associated to distinct families of structure group constants 

C:, = eaMN''' + S',A,-S:Ab (1) 

where Eabd is the usual totally antisymmetric unit tensor, 6^ the Kronecker symbol, N"-^ 
is the contravariant component of an order 2 symmetric tensor, and the vector A must 
follows 

N'^'^Ab = NabA'' = 0. (2) 

Without loss generahty, one can write A^"** = diag(ni, ^2, ria) and A^ — [a, 0, 0] provided 
that ani = 0. Distinct homogeneous and anisotropic spaces in 3D can then be classified 
in the following table: 
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Bill 


is not an eigenvalue of N 
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Bix 
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1 


-1 





Byiii 



Note that each case is degenerate, for example Ui — —1, n2 — 1, — —1 and a — 
is in the equivalence class of Bix which contains all the possibilities with a positive 
signature and out of the spectrum of N. This is the well known Bianchi classification 
(see [1],[2]). 



1.2. Scalar-tensor theory of gravitation 

In scalar-tensor theories of gravity, the dynamics of the Universe contains a new 
scalar degree of freedom that couples explicitly to the energy content of the Universe 
[3, 4, 5, 6, 7]. In units of c = 1, the action generically writes, in the so-called Einstein 
frame: 

+ S„,{^IJrn:Q^{^)9nu): (3) 
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G being a bare gravitational constant, ip the scalar field, U{ip) its self-interaction term 
and 0(v?) its coupling to matter. The functional S„i{ijyni- ^'^{v)9^lu) stands for the action 
of any field ipm that contributes to the energy content of the Universe. It expresses 
the fact that all these fields couple universally to a conformal metric g^i, = 0^((/9)(?^j,, 
then implying that the weak equivalence principle (local universality of free fall for 
non-gravitationally bound objects) holds in this class of theories. The metric g^,, 
defines the Dicke- Jordan frame, in which standard rods and clocks can be used to 
make measurements (since in this frame, the matter part of the action acquires its 
standard form). Despite the conformal relation, these two frames have a different 
status: in the Dicke- Jordan frame, where the gravitational degrees of freedom are 
mixed, the Lagrangian for the matter fields does not contain explicitly the new scalar 
field: the non gravitational physics has then its standard form. In the Einstein frame, 
the scalar degree of freedom explicitly couples to the matter fields, then leading for 
example to the variation of the inertial masses of point-like particles. Of course, the two 
frames describe the same physical world. Nevertheless, the usual interpretation of the 
observable quantities is profoundly modified in the Einstein frame, whereas it holds in 
the Dicke- Jordan frame, where the rods and clocks made with matter are not affected 
by the presence of the scalar field. That is why one usually refers to the Dicke-Jordan 
frame as the observable one. However, the dynamics of the fields is generally more 
easily described in the Einstein frame, so that in this work, since we are interested in 
the integrability of the models rather than in there physical content, the analysis will 
be done in the Einstein frame. 

Varying the Einstein frame action (3) with respect to the fields yields the equations: 

i?M--^%^. = xi;.+ij. (4) 

n^ = --ui^)T + ^ (5) 

V.r; = u;{ip)TV^ip (6) 

where x = SttG, T is the trace of the energy-momentum tensor of matter fields T^,,, 
and = 2if, (g'^'^LPaVls) - '2U {^)gf,u is the energy-momentum tensor of the 

scalar field. Moreover, we have defined the coupling uJ {(f) = It is important to 

note that these equations reduce to those of General Relativity in presence of a scalar 
field iff a; (99) = 0. 

1.3. Dynamical equations of scalar-tensor theory in homogeneous spaces 

If one denotes by t the physical time and by r a conformal time such that dt = N (r) rfr, 
where N (r) is generally called "lapse function", the line element of the 3-1-1 physical 
space reads: 

ds^ = g.^dx'dx^ -dt^ = -f (r) uj'uj^ - (r) dr^ (7) 

As explained in [8], one can find an invariant basis of differential forms {u;^, u;^, cu^} in 
each Bianchi space such that 7 is diagonal; we note hereafter 7 (.) = diag (a^ (.)) . 
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Taking V (r) — {aia2a^y^^ as a lapse function, and introducing Ai — In (ctj), the 
scalar-tensor dynamical equations in homogeneous and anisotropic Universe read: 



[p + p^-iP + P^)] 
XV [p + p^-{P + P^)] 



A'l + njaf - {n2a2 - 7130:3) 
A2 + n^al — (nstts — nicti)^ 
^3 + n|Q;| — (nitti — 77,20:2)^ 
-2(^"-4(lnl^)V'-2^l^' 



(8) 



In these equations we have written ' for d/dr; in addition the Universe is filled by 
a perfect fluid with pressure P and energy density p; finally we have noted: 

: b'V2 + f/(^)] /X 



P^ := [^'y2-U{^)]/x 



(9) 



Taking U = 0, and reorganizing the first dynamical equation using the three others, 
the system (8) becomes: 



= 



XV [p + p^-{P + P^)] 
XV' [p + p,^-{P + P^)] 

XV [p + p^-iP + P^)] 
xV'u {^) [p - 3P] 



Ec + Ep- AxpV - 
A'[ + n\a\ - {n2a2 



^XP^V'^ 



A2 + H2a| — (na^a — niaiY 
A3 + nlal — {riiai — n2Q!2)' 
-2(p" - 4 (In V)' cp' 



(10) 



where 



E 

1=1 



2„2Ai 



n, e 



:= a; A'2 + A;A^ + A'3A'2 and Ej,:=Y^ 

The energy-momentum conservation, that in scalar-tensor theory is: 

W^T^.^uiip)T^fsg^^d,ip (11) 
coupled to the hypothesis of a barotropic fluid: 

p = (r-i)p (12) 

allows us to obtain a relation between p, V and (f which is: 

p = p^y-^e^-^^ with po e R+ \ {0} (13) 

One should note that a fluid of radiation (F = 4/3) doesn't couple directly to the scalar 
fleld. Assuming a power law dependence of </?' in V, more precisely: 

= v'oV^ (14) 

and making use of (13) one can solve the last equation of system (8), and obtain an 
explicit dependence on V for 6(<^): 

4(A + 2) 



p^Q^-^^ = ^^2(A-l)+r ^j^g^g ^2 



2X 



2(1- A) -r 



(15) 
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Assuming > (which corresponds to non exotic matter) and < F < 2 (which 
corresponds to the largest class of barotropic fluid), relation (15) holds under the 
condition: 

-2 < A < ^ (16) 

The assumption (14) is a constraint on the entire dynamical system rather than on 
the scalar-tensor theory itself. Indeed, any choice of coupling function u!{ip) can be done, 
but then the resulting behavior of V{t) is completely fixed by the last equation of system 
(10). Conversely, imposing a behavior for the volume V{t) determines the corresponding 
coupling function. For example, imposing a Brans-Dicke theory, i.e. Lj((f) = ujq = cste 
results in V{t) oc r^/*^^~^); on the contrary, a volume evolving asVoc e^'^ leads to 

2(A-i)+r 

e((/7) oc (p (4-3r)A ^ (17) 

that is, u{(p) oc ^ and ip oc e^^'^. To sum up, the only constraint imposed by assumption 
(14) is on the couple {V{t),u!{(p)), through the relation: 

A = ,i(^v- . (18) 

2(1 - A) - F 

that implies that, at any time, the ratio of the densities of the barotropic fluid and 
of the scalar field scales with the inverse of the square of the volume. This is a 
sufficient constraint to make the Kovalewski formalism tractable in scalar-tensor gravity. 
Replacing (15) into (13), and considering (14) after (9), the dynamical system associated 
with the homogeneous Universe in scalar-tensor theory is: 

E, + E^ = 2a2\/2A + 2v9;V2^+2 

2[2-F]a2\/2^ = A'; + nlo^~{n2a2-nsasf 

2[2-F]aV2^ = A'^ + nlal - {n^a^ - n,a,f ^ ^ 

2[2-F]a2\/2A _ A'^ + njaj - (niai - n2a2f 

It is important to note that whereas this system seems independent on F, it actually 
depends on it through the parameter A that fully characterizes our solution. In fact, 
for each coupling function, their exists a non-ambiguous link between F and A. For 
example in the case of a radiation fiuid when F = 4/3, the LHS of the last equation in 
system (10) vanishes and one can find that cp' oc V~'^ and then A = —2. 

A direct inspection of the first equation of system (19) shows that we can predict 
quahtatively the behavior of the case: A = 0. Indeed, in this case, one of the two terms 
of Ec + Ep reduces to a positive constant and the other one tends to as K tends to 0. 
It is well known (e.g. [10]) that such a case breaks the Kasner cycle (e.g. [11]) for -Byill 
and BiXj and then suppresses the chaotic behaviour toward the t = singularity for 
these models. The case A = — 1 seems to be similar, but the other term now diverges 
and then this simple analysis cannot be done. 
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1.4- Hamiltonian formalism 

The quantity called Ec — A[A2 + A[A'^ + A'2A2 is a quadratic form of >li=i,2,3 derivatives. 
Then, one can diagonalize it using a linear change of variables: 

q, = {A,-A2)/V2 

q2 = (^1 + A2- 2^3) /\/6 (20) 
53 = 2 {Ai + A2 + As) /VQ 

Introducing the associated conformal time derivatives Pi=i,2,3 := ^'^^ 23' ^^^^ equation 
of the system (19) becomes: 

H := I (p^ -p? -p^)+e^«3^(^^^^^)_2a2e^A,3_2<^;2^^(A+i),3 = o(21) 
the so-called potential ^ is defined by: 

e (gi, 52) = -n?e(^^^+^^0 - n^e(^^-^^0 - n^e"^^^ (22) 

+ 2n,n2e'^^^ + 2nin3e(^^-^'^0 + 2n2n3e-(^^^+^''^) 

In terms of {qi,Pi) variables, the dynamical system (19) is quasi-Hamiltonian (e.g. 
[12] and later [13], [14], [15]) 

/ _ rfgi,2 _ dH I _ dpi,2 _ dH 
"1>2 dT dpi 2 "1,2 dr dq\ 2 

(23) 

"3 dT dps "3 dr dqz 

The minus that breaks the strict Hamiltonian symmetry comes from the minus signature 
of the quadratic form associated with Ec- 



2. Integrability of homogeneous Universes 

2.1. The case of Bianchi Universes 

The following work has been initiated by Melnikov's team (see e.g. [16] and [17] 
and references therein). In the special case of -Bjx, [18] consists in an application; a 
generahzation in the context of the whole class A (i.e. a = 0) was tried by [19]. However, 
a lot of imprecisions in this last work need this new reformulation and extension to 
scalar-tensor theory. 



2.1.1. Bianchi universes as generalized Toda systems Introducing the following vectors: 



ai 
a4 
ar 



= [0, v^/3, v^/3] a2 := [^2/2, v^/3] a3 

= [V2, V6/3, V6/3] as := [- V2, V6/3, V6/3] ag 
= [0, 0, V6A/2] ag 



the 3-forms: 



Vx,y e 



(x, y) := +xiyi + X2y2 + x^y^ 
(x, y) := -xiyi - X2y2 + x^ys 



= [-y2/2,-v^/6,v^/3] 

= [0, -2v^/3, V6/3] (24) 

= [0,0,V6(A + l)/2] 



(25) 
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and the constants: 

ki :— 2nin2 ^2 := 2nin3 k^ := 277,2713 



k 



4 •- 



-77,; 



kf. 



-nl 



ke, :— —n% 



(26) 



k7 = -2a^ k^ = -2(f'^ 
it is clear that the Hamiltonian (21) reads: 



1 ^ 



(27) 



i=l 



which is a classical form of a generalized Toda dynamical system. 

Following Melnikov [16], we change {q, p} variables to {u, v} ones, such that: 

u e Mj=i,...,8 (a», p) 
V e R^, Vi=i^,„^8 := exp (a^, q) 

Through this change, the number of degrees of freedom jumps from 6 in {q, p} to 16 in 
terms of {u, v}. Still writing ' for the derivatives with respect to the conformal time r, 
the dynamical equations then read: 



(28) 



UiVi 



with rriij :— —kj {sii, aj 



(29) 



TUijVj 



i=i 

This new formulation is now polynomial. Using the appendix notations, one can directly 
prove that the system (29) is autosimilar with any non vanishing index and weight g 
such that: 

gi = ... = gs = l and gg = ... = g^Q = 2 (30) 

g is unique provided that (^.3) is fulfilled. A particular autosimilar solution of (29) is 
then 

[uv]T = crs= [Xit-\...,X8t-\iiit-^,...,fi8t-^f (31) 

provided that the constant non vanishing vector c = [A, /i] = [Ai, As, //i, //g] is a 
solution of the algebraic system of equations: 



Vi = 1, 



-A,; 



(32) 



Ai — — 2//j 

Hence, to any non vanishing solution of this last system corresponds a set of 16 
Kovalewski exponents that allows to write the solution of the system (23) (cf Appendix 
A). A necessary condition for the system to be integrablc is that all its Kovalewski 
exponents be rational. So, the rest of the paper will be devoted to the analysis of these 
exponents in order to study the integrability of different types of homogeneous and 
anisotropic Universes. 
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2.1.2. Solutions of the algebraic system In the more general case (that is scalar-tensor 
theory in presence of matter barotropic fluids), the algebraic system (32) makes use of 
the 8x8 matrix 








— 2nin3 


-2n3n2 








2ni 


20,2 A 


2v^;^Ai 




— 2nin2 





-2n3n2 





2nl 





2a2A 


2¥^;^Ai 




—2n\n2 


—2n\nz 





2n? 








2a2A 


2<^:,^Ai 










— 4n3n2 


-2nf 


2n2 


2n| 


2a2A 




M := 





— 4nin3 





2nl 


-2n^ 


2n| 


2a2A 


2v^'o'Ai 




— 4nin2 


— 4nira3 





2nl 


2ni 


-2ni 


2a2A 


2v^'o'Ai 




— 2nin2A 


— 2nin3A 


-2n3n2A 


nfA 


niA 


niA 


Sa^A^ 


3(^;2AiA 




— 2nin2Ai 




— 2n3n2Ai 


nfAi 


n|Ai 


nlAi 


Sa^AAi 





V(i, j, k) eEpXEpX Ep Q = , C 



,i,j,k 



m 
m 



ITii 



ma 
rriji 
rriki 



m. 



rriik 
rrijk 
rrikj rukk 



m 



33 



m 



33 



(34) 



and the determinants: 



d2i = 



rrii 



m 



'33 



m 



m 



'33 



rriik 
rrijk 



rrikj rukk 



d2j 




ma 




rriji 


ma 


2 


rriik 


rriji 


2 


rrijk 


rriki 


2 


rukk 



d: 



3fc 



ma 
rriji 
rriki 



rrii 



m 



'33 



rrikj 



2 
2 
2 



Solutions of System (32) arc then classified into three classes: 

(i) Type 1 solutions (Tl): 3\ i E Ep such that Hi ^ and Vj E Ep \ {i} , fij 
A, = — 2 and 



(35) 



(36) 



: then 



• if Ci = : there is no solution. 

• if Ci 7^ : /ij = 2/Ci and Vj e Ep\ {i} , \j = -2mji/(i 

(ii) Type 2 solutions (T2): 3! {i,j) e Ep x {Ep \ {i}) such that {iii^i^j} 7^ {0,0} and 
Wk E Ep \ {i,j} , fik = '■ then Aj = \j = —2 and 

• if (ij = : there is no solution. 

• if (ij 7^ : fii = d2i/Cij and /^j = d2j/Ci,j moreover \/k e Ep \ {i,j}, 

Afc = irnkid2i + mkjd2j) /Qj 

(iii) Type 3 solutions (T3): 3! {i,j,k) e Ep x{Ep\{i}) x {Ep\{i,j}) such that 
{/Ii, iij, Ilk} 7^ {0, 0, 0} and Wl G Ep \ {i,j, k} , /ii = : then Aj = Xj = A^ = — 2 and 



(33) 



where Ai := A + 1. One can straightforwardly check that M is of rank 3. 

The vector is always a solution of the system (32), but it is not relevant for the 
quest of Kovalewski exponents that needs non trivial solutions. A systematic study of 
the solutions of the system (32) is possible : for p = 6, 7 or 8 let Ep = {1, 2, 3, • • • 
we consider the minor determinants: 
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if Ci,j,k 



= : there is no solution 



d3j/Ci,j,k and /i^ = d3k/Ci,j,k moreover 



VI e Ep\ {i,j, k}, Xi = {miidsi + rriijdsj + mikd^k) /Ci,j,k 



2.1.3. Kovalewski exponents for Bianchi Universes In what follows, we will examine 
the Kovalewski exponents for Bianchi Universes in 4 different cases that are included in 
the formalism presented above: scalar-tensor gravity with or without matter fluids, and 
General Relativity with or without matter fluids. We proceed in three steps. 

(i) Select a Universe. It corresponds to choosing a set of nj=i^2,3 in table 1.1. This step 
determines the theory of gravitation and the matter content of interest, and this 
choice is made by considering different forms for the matrix M defined in (i): 

• Barotropic matter filled Universe in scalar-tensor theory of gravitation: p — 8. 
This is the most general case that was presented in detail in the preceding 

subsection. 

• Empty Universe in scalar-tensor theory of gravitation: p = 7. The associated 
matrix is M but without the 7th line and the 7th column: 








— 2rain3 


-2n3n2 








2nl 






— 2nin2 





— 2n3n2 





2nl 





2<^'.^Ai 




— 2nin2 


— 2nin3 





2nl 








2^'o'Ai 


MsTV '■ — 








-4n3n2 


-2nl 


2n\ 


2nl 


29^'o'Ai 







—Auius 





2n\ 


-2n\ 


2n\ 


2v^'.2Ai 




— 4nin2 


— 4nin3 





2n\ 


2nl 


-2nl 


2<^;^Ai 




— 2nin2Ai 


— 2nin3Ai 


— 2n3n2Ai 


n?Ai 




niAi 


3<^;^Af J 



Barotropic matter filled Universe in Einstein General Relativity: p — 1 . The 
associated matrix is M but with the 8th line and 8th column removed, and 
the 7th line and 7th column adapted. Noting 7 for 2 — F we have: 








— 2nin3 


-277,3712 








2nl 






— 2nin2 





— 2n3n2 





2ni 





'^PoXl 




— 2nin2 


— 2nin3 





2nl 








'^PoXl 


Mqrm '■— 








-4n3n2 


-2n\ 


2nl 


2ni 


'^PoXl 







— 4rain3 





2n\ 


-2n\ 


2nl 


'^PoXl 




— 4nin2 


— 4riira3 





2n\ 


2n\ 


-2nl 


"^PoXl 




. — nin27 


-nin37 


—7137127 


n\ll2 


nh/2 


nh/2 


^poxiV^ 



M, 



GRV ■'- 



n 



Empty Universe in Einstein General Relativity: p — Q>. The associated matrix 
is M but with the 7th and 8th lines and the 7th and 8th columns removed: 


— 2711772 

— 27li7l2 




-471x712 

We then have to apply the following algorithm to the appropriate matrix. 

Determine all the non vanishing minor determinants C, which could be extracted 
from the considered matrix. 



277-1713 


-27l3n2 








2nl 





-2nzn2 





2nl 





27Xi7l3 





2nl 











—4713712 


-2nl 


2ni 


2ni 







2nl 


-2ni 


2"-3 


4nin3 





2nl 


2nl 


-277 
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(iii) For each C, compute the associated set of 2p Kovalewski exponents. In our 
polynomial indicated in appendix A this set is the set of eigenvalues of 

the matrix: 

1 mil ?7ipi 

■•. ■•. ; ; ; 



K 



M 



••• 

/^i 

■-. 







■-. : 

1 mip 

■■■ Ai + 2 

: 







VP 











Ap + 2 



(37) 



■-. 

iJip 

The particular form of the matrix K and the weak rank of the matrix M allows 
an explicit calculation of the eigenvalues in all the cases p — 6, 7, 8 and for all the 
Bianchi Universes. There exists several hundreds of non vanishing minor determinants 

in all the cases considered. They correspond to 92 distinct sets of Kovalewski exponents 
with length 2p = 12, 14 or 16. All the details arc presented in appendix B. In all cases, 
calculations are explicit, and in order to organize it we use formal calculus computational 
tools. 



2.1.^. Exponents analysis and conclusions The full set of Kovalewski exponents is 
presented in appendix B. The analysis of the sets of exponents results in four classes 
for Bianchi models : Class I contains uniquely models. Class II contains Sjj and 
Sjy, Class III is composed by -Bjjj, -Byj and -Byjj , and finally Class IV contains 
Syjjj and the famous -Bjx- Each class is characterized by the same set of Kovalewski 
exponents, therefore corresponding models have the same characteristic dynamics. Let 
us review for each class the algebraic properties of Kovalewski exponents : 

• Class I — For all the cases we studied, all the Kovalewski exponents belonging to 
this class, are rational provided that the barotropic index V and/or the scalar-tensor 
parameter A are rational. This restriction corresponds to physical casesf. 

• Class II — Kovalewski exponents of models belonging to this class fall into four 

subcases : 

— Empty Universe in General Relativity (EUGR): exponents belonging to this 

class are integers. 

— Barotropic Matter filled Universe in General Relativity (BMUGR): due to A*^, 
all exponents are rational iff F e Q n [0, Fq] with Fq = (11 - V73)/3 ^ 0.81. 

I As a matter of fact, physical power laws or barotropic index must be rational in order to be full of 
physical meaning 
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— Empty Universe in Scalar- Tensor Theory (EUSTT): due to all exponents 
are rational iff A e Q n [Aq, 1] with Aq = (-11 + ^)/6 ^ -0.40. 

— Barotropic Matter filled Universe in Scalar- Tensor Theory (BMUSTT): Due 
to A± and D±, all exponents are rational iff A e Q n [Ai, 1] with Ai = 
(-5 + v^)/6 ^ 0.59. 

• Class III — As previously, four classes of sets of Kovalewski exponents can be 
identified 

— EUGR : exponents belonging to this class are integers. 

— BMUGR : due to A*^, B*^ and C^, all exponents are rational iff T e Qn [0, Tq]. 
This case is then equivalent to the corresponding case of class II models. Let 
us remark the special value V = 2/3 for which {A\, A*_, B*^, B*_, C% C*_) = 
(1, 0, 1, 0, 1, 0). 

— EUSTT : due to A±,B± and a± all exponents are rational iff A e Q fl [Aq, 1]. 
This case is then equivalent to the corresponding case of class II models. Let us 
remark the special value A = —1/3 for which {A^, A_, B^, B_, a+, a_) = 
(1, 0, 1, 0, 1, 0). 

— BMUSTT : due to A±,B±, C±, D±, a± and b±, all exponents are rational iff 
A e Q n [As, 1] with As = 16/25 = 0.64. 

• Class IV — all models of this class contain at least two conjugated complex 
Kovalewski exponents. 

/^From such an analysis, taking into account Yoshida's theorems (see [21, 22]) whose 
context is detailed in appendix A, we can conclude that: 

• Empty Universes whose metrics correspond to classes I, II or III defined below are 
generically associated to algebraically integrable dynamics. 

• Class IV Universes are always associated to non algebraically integrable dynamics. 
This result holds for General Relativity and/or scalar-tensor theory we are 
interested in, the Universe being empty and/or filled of barotropic matter. 

• BMUGR of class II and III arc associated to non integrable dynamics if the 
barotropic index ranges in the interval [ro,2] with Fq = (11 — V73)/3 ~ 0.81. 

• EUSTT of class II and III are associated to non integrable dynamics if the power 
law of the scalar-tensor modelization (14) ranges in the interval [— 2, Aq] with 
Ao = (-11 + V73)/6 ^ -0.40. 

• BMUSTT of class II are associated to non integrable dynamics if the power law 
of the scalar tensor modehzation (see 14) ranges in the interval [—2, Ai] with 
Ai = (-5 + V73)/6 fa -0.40. 

• BMUSTT of class III are associated to non integrable dynamics if the power 
law of the scalar tensor modelization (14) ranges in the interval [—2, A2] with 
A2 = 16/25 = 0.64. 

The integrability of homogeneous Universes depends on multiple factors. 
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It is well known, since the pioneering works by Elskcns and Henneaux [23], that the 
number of dimensions D of the spatial sections of the Universe is one of this factors. It 
influences the form of the potential ^ (q) in relation (22) or/and increases the number 
of components of the vectors {aj} in (24). The net result in the simplified context of 
generahzed Kasner metric is that chaos disappears — and then system becomes always 
integrable, in a classical context — when D > 11. 

Modification of the gravitation theory by introducing a scalar-tensor component 
seems here to have a different contribution. As a matter of fact, in the case we 
consider, the potential ^ (q) is not modified, but a volume dependent Toda potential 
is added in the Hamiltonian. In a general way, this produces a slowing down for 
the cushions' velocity in the billiard representation (see [15] and references therein). 
One can then reasonably suppose that this feature does not change drastically the 
dynamical properties of homogeneous Universes inherited from General Relativity. The 
integrability analysis based on Kovalewski exponents we have done confirms such a 
conjecture. 

Appendix A : Integrability of autosimilar differential systems 

Let f be a continuous function from M" to M" and x = [xi{t), Xn{t)]^ G M" such that 
^-f(x) (A.1) 

If there exists a positive real A and a rational vector g :— [gi,...,gn] such that the 
transformation 

^ ^ (A.2) 

Vi = 1, n X^^Xi 

leaves the system (^.1) invariant, this system is called autosimilar with index A and 
weight g. 

When it exists g, it is unique if 

df - fx) 

A (x) e M„ (M) , A,j = - Sijfi (x) (A.3) 

is invertible for almost all x eR". 

Autosimilar systems always admit at least one autosimilar particular solution 

ias= [ci{t-to)-'\...,Cn{t-toy'"f (A.4) 

where c := [ci, .... c„]^is solution of the algebraic equation: 
/i (c) = -gici 

: (A.5) 

fn (*-') 9n^n 

Linearization of the system (A.l) around the solution Xa^ gives 

^ = Di {^) (SCas) z (A.6) 
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where Df (x) (x^s) is the Jacobian matrix of f (x) evaluated at x = x^s- A theorem by 
Fuchs (see [20] §5.6 and 5.7) then shows that the general solution of (A. 6) is autosimilar 
and reads: 

z = [^1 {t - tX''' , K {t - t^r-'-f (A.7) 

The quantities {pi, p2, Pn) are called Kovalewski exponents. More practically, one can 
compute them because there are also eigenvalues of the matrix: 

K := Df (x) (c) + diag (g) . (A.8) 

As shown by Comte (see [20] §5.6 and 5.7), they are of great importance for the 
study of integrability of the original non linear system {A.l). As a matter of fact each 
component Xi=i_..^„ of a solution x of this system can be written as: 

Xj^ Ci (t to) 



hoo 



Xi{t)^J2e''x\''\t) with 



k=0 



(2) _ „ J, /j- 7- \9p+Pq-gq 



3 1 < p,q < n, xy = Cpkq {t - to) 
etc • • • 



hence 

X, (t) oc{t- toT s [{t - toY' , {t - toY-] (A.IO) 

where 5* [ ] is a multiple series. This allows to understand a theorem by Yoshida (see 
[21, 22]) : a necessary condition for a differential system to be integrable is that all its 
Kovalewski exponents are rational. If there exists at least one exponent irrational or 
complex the corresponding differential system is not algebraically integrable. 

Appendix B : Details of Kovalewski exponents 

When only non exotic barotropic matter ( < F < 2) is considered, we have noted 
7 = (3 F - 2)/(F - 2). For any X = 3, 5,7 and 15 we note 2K^^ = 1 ± i^fK and 



2A1 = 1 ± ./ le-fr+Br^ 2S1 = 1 ± y/b^^ 



2Ci = 1 ± ^J^E^ 2Dl = l± y^ i8-4ir+24r^ 
£;^:={-l,2, l(x6),f (x6)} £;„:={-l,4(x3),2(x3),l(x5)} 

Moreover 

^5= [{-1,1 (x6) ,2 , 2V3(x6)} 

(5i = (3A + 1) (A + 1)-^ and ^2 = (3 A - 2) A"^ 



/3(A- A+) (A- At) , ^ 

2A±^l±\j^ ^-^p^- ^ 6A^ = -ll±v^ 
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2D± = 1 ± 



'3 (A - A+) (A - A^) 



6A^ = -5 ± \/73 



2B± = 1 ± 



'25A + 9 
|A + 1| 



2C+ = 1 ± 



/25A- 16 



2£;± = 1 ± 



/48A2 + 41A + 9 



lA + ll 



48A^ + 41A + 9 > 



2F± = 1 ± 



/48A2 - 55A + 16 



48A^ - 55A + 16 > 



2a± = 1 ± V5 + 12 A 26± = 1 ± V-7 + 12 A 



Using the notations defined above we have computed all the Kowalewski exponents 
for all Bianchi Universes, in the case of empty or barotropic filled Universes and for 
General RcIativity(GR) and Scalar-Tensor (ST) theory of gravitation: 





RG 


Empty (N=6) 


Barotropic Matter (N=7) 


Tl 







1'2 








'1'3 








ST 


Empty (N=7) 


Barotropic Matter (N=8) 


Tl 


Eg 


-1, , l(x7), 2 , 252/3(x6), -2/A 
-1. l(x7).2.2(A^l)-^ 2(5i/3(x6) 


12 








1'3 









Bii and Diy 


GR. 


Eiupty (X=6) 


Barotropic Matter (X=7) 


Tl 


Eo 


-l,l(x6), 2(x3), 4(x3),3-r/2 


T2 





-l,l(x5),2,7/2(x2),7(x3),A| 


'1'3 








ST 


Empty (N=7) 


Barotropic Matter (N=8) 


Tl 


-1, l{x6), 2(x3), 4(x3),A + 3 

Eg 


-1, l{x7),2, 2(A + l)-\2<5i/3(x6) 

-1, l(x7), 2,252/3{x6), -2/A 

-1, l(x7), 2(x3), 4(x3), A + 3, A + 2 


T2 


-1, l(x5), 2, <5i/2(x2), <5i(x3), A± 


-1, l(x6), 2, fc/2(x2), 52(x3) , -2/A, D± 

-1, l(x6), 2, 5i/2(x2), <5i(x3), 2(A + l)-\ A± 


'1'3 
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Bill , /„,,,, Bvii^_^ 


GR 


Empty (N=6) 


Barotropic Matter(N=7) 


Tl 


Eo 


-1 , l(x6), 2(x3), 4(x3), 3 - r/2 


T2 





-l,l(x5),2,7/2(x2),7(x3),A5^ 
-1, 0(x2), l(x5), 2, 27, 7(x2), 


T3 





-1, 0, l(x4), 2, 27, 7(x2), SJ, 
-1, 0(X2), 1(X5), 2, 27, 7{x2), 


ST 


Empty {JN=7) 


Barotropic Matter (N=8) 


Tl 


— 1, l(x6), 2(x3), 4(x3),A + 3 


-1, l(x7),2, 25i/3(x6), 2{A + 1)"-' 

-1, l(x7), 2,2<52/3{x6), -2A-i 

-1, l{x7), 2(x3), 4(x3),A + 3, A + 2 


T2 


-1, l(x5), 2 , 5i/2 (x2), 5i(x3), A± 
-1, 0(x2),l(x5), 2,2<5i, <5i(x2), a± 


-1, 0(x2), l(x6), 2, 2<5i, 5i(x2), 2{A + l)~\a± 
-1, l(x6), 2, 5i/2 (x2),5i(x3), 2(A + 1)~\ yl± 
-1, 0(x2), l(x6),2, 252, 52{x2),-2A-i, b± 
-1, l(x6), 2, <52/2{x2), 52{x3),-2A-i, D± 


T3 


-1, 0, l(x4), 2, 25i, <Si(x2), a±, B± 
-1, 0(x2), l(x5),2, 25i, <5i(x2), a± 


-1, 0, l{x5) ,2, 26i, 5i(x2), 2{A + a±, B± 
-1, 0(x2), l(x6), 2, 2(5i, <5i(x2), 2(A + 1)"\ a± 
-1, 0(x2), l(x6), 2,252,<52(x2), -2A-1, fe± 
-1,0, l(x5), 2,52(x2), 252,-2A-i,6±, C± 



Bviii and Six 



GR 


Empty (N=6) 


Barotropic Matter {N=7) 


Tl 


Eo 


-1, l(x6), 2(x3), 4(x3),3-r/2 


T2 


-l,-2(x3),7±% 2(x2), l(x4) 
-1, 2, 15±% l(x4), 0(x4) 


-1, 1 (x5) , 2, 7/2 (x2) , 7 (x3) , yl5^ 
-1, 0(x2), l(x5), 2, 27, 7{x2), 
-2(x3),-l, l{x5), 2(x2),7±*,2r- 2 
-1, 0(x4), l(x5), 2, 15±*,-H-3r/2 


T3 


-1,2, 15±'(x2), l(x3), 0(x3) 
-1,2, 3±'(x2), l(x3), 0(x3) 
-1, 2, 15±', l{x4), 0{x4) 
-l(x2),2(x2),7±S l(x3), 0(x3) 


-1, 0(x3), l(x4), 2,3±'(x2), -1 + 3r/2 

-1, l(x4), 2, -27, 27(x3),C;!;^,D5^ 

-1, 0(x3), l(x4), 2, 15±'(x2),-l + 3r/2 

-1, 0(x4), l(x5), 2, 15±',-l + 3r/2 

-1, 0, l(x4), 2, 27, 7(x2), BJ, 

-l(x2), 0(x3), l(x4), 2(x2), 7±%-l + 3r/2 

-1, 0(x2), l(x5), 2, 27, 7(x2), B% 


ST 


Empty (N=7) 


Barotropic Matter (N=8) 


Tl 


Es 

-1, l(x6), 2(x3), 4(x3),A + 3 


-1, l(x7), 2,252/3(x6), -2A-i 

-1, l(x7), 2(x3), 4(x3), A + 3, A + 2 

-1, I(x7),2,25i/3(x6), 2(A + l)~i 



T2 



-1, 0(x2), l(x5), 2, <5i{x2), 25i, a± 
-2(x3),-l, l{x5), 2{x2), 7±*,-4A-2 
-1, 0(x4), l(x5), 2, 15±*,-3A - 1 
-1, l(x5), 2,<5i/2(x2), <5i(x3),A± 



-1, 0(x2), l{x6), 2,2<52, 52(x2), -2A-\ b± 
-1, 0(x2), l(x6), 2, 2<5i, 5i(x2),2(A + l)"\ a± 
-1, l(x6), 2, 5i/2(x2),<Si(x3), 2 (A + 1)"^ ,A± 
-2(x3),-l, l(x6), 2{x2), 7±',-4A- 2, -4A + 2 
-1, 0(x4), l(x6), 2, 15^',-3A- 1, -3A + 2 
-1, l(x5), 2, -25i, 26i(x3),2(A + l)-\ B±, E± 
-1, l(x5), 2,-252, 252(x3),-2A-i, C±, F± 
-1, 0(x3), l(x5), 2, 3±'(x2),-3A- l,-3A + 2 
-1, 0(x3), l{x5), 2, 15±'(x2),-3A- 1, -3A + 2 
-l(x2), 0(x3), l(x5), 2(x2), 7±% -3 A - 1, -3 A + : 
-1, 0, l(x5),2, 25i, <5i(x2), 2(A + 1)"\ a±, B± 
'1,0, l(x5), 2, 2(52, 52{x2), -2A-1, 6±, C± 
-1, 0(x2), l(x6), 2,2(52, (52(x2),-2A-i, b± 
-1, 0{x2), l(x6), 2, 2(5i, (5i(x2), 2(A + 1)~\ a± 
-1, 0(x4), l(x6), 2, 15±%-3A- 1, -3 A + 2 



T3 



-1, 0(x2), l(x5), 2, 2<5i, 5i (x2), a± 
-1, 0(x4), l{x5), 2, 15±',-3A- 1 
-1, 0(x3), l{x4), 2,15±^(x2), -3A-1 
-1, 0(x3), l(x4), 2, 3±'(x2),-3A- 1 
-1, 0, l(x4), 2, 2(5i,(5i(x2),a±. B± 
-l(x2), 0(x3), I(x4),2(x2), 7±',-3A-l 
-1, l(x4), 2,-25i, 25i(x3), B±,E± 
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